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ABSTRACT 



o 
o 



^ ' We present a model for the non-thermal radio emission from presumably single O stars, in terms of synchrotron emission from relativistic 
f — ' electrons accelerated in wind-embedded shocks. These shocks are associated with an unstable, chaotic wind. The main improvement with 
CN| , respect to earlier models is the inclusion of the radial dependence of the shock velocity jump and compression ratio, based on one-dimensional 
' time-dependent hydrodynamical simulations. The decrease of the velocity jump and the compression ratio as a function of radius produces 
I a rapidly decreasing synchrotron emissivity. This effectively prohibits the models from reproducing the spectral shape of the observed non- 
thermal radio emission. We investigate a number of "escape routes" by which the hydrodynamical predictions might be reconciled with the 
radio observations. We find that the observed spectral shape can be reproduced by a slower decline of the compression ratio and the velocity 
, jump, by the re-acceleration of electrons in many shocks or by adopting a lower mass-loss rate. However, none of these escape routes are 
• physically plausible. In particular, re-acceleration by feeding an electron distribution through a number of shocks, is in contradiction with 
Oh' current hydrodynamical simulations. These hydrodynamical simulations have their limitations, most notably the use of one-dimensionality. At 
^ , present, it is not feasible to perform two-dimensional simulations of the wind out to the distances required for synchrotron-emission models. 
5^ ' Based on the current hydrodynamic models, we suspect that the observed non-thermal radio emission from O stars cannot be explained by 
' wind-embedded shocks associated with the instability of the line-driving mechanism. The most likely alternative mechanism is synchrotron 
Ci • emission from colliding winds. That would imply that all O stars with non-thermal radio emission should be members of binary or multiple 
^ ' systems. 
• I— I 

' Key words, stars: early-type - stars: winds, outflows - radiation mechanisms: non-thermal - hydrodynamics - instabilities - binaries: general 

^' 

_ _ I 

1 . Introduction by a factor of ten at 6 cm by Bieging et aL ll989> . It also has a 

very different spectral shape, with a negative spectral index, 
i.e. a flux that decreases as a function of frequency. Such a 
spectrum is characteristic of synchrotron emission from rela- 
tivistic electrons gyrating in a magnetic field. Electrons can at- 
winds, they are all expected to have thermal emission, although ^^^^ relativistic speeds through first order Fermi acceleration in 
It IS intrinsicaUy faint and can only be detected for the bright- ^j^^^j.^ (g^jj ^973 Drury 1983;. Although shock-acceleration 
est and closest stars. A significant fraction of hot stars also j^^^ ^^^^ ^^^^ frequendy used to explain synchrotron emis- 
emits non-thermal radiation, which is generaUy attributed to ^j^^ f^^^ ^^^^^^ ^^^^^ mechanisms such as magnetic recon- 
synchrotron emission from relativistic electrons. „g^jion are also capable of accelerating particles to relativistic 

Wright & Barlow ( I1975> developed an elegant formalism speeds (e.g. Litvinenko lMBt . 
for the thermal radio emission from hot stars. Assuming a uni- 
form, isothermal, spherically symmetric outflow, they showed 

that the expected flux is of the form Fy oc v", where the spectral A fundamental question regarding non-thermal radio emis- 
index a - +0.6. The same formalism aUows a star's mass-loss sion from O stars is its correlation with binarity. If the non- 
rate to be derived from its thermal radio flux. thermal emitter is a member of a binary system, the shocks 
Non-thermal emission from hot stars is generally much needed to accelerate the electrons can be provided by colliding 
stronger than thermal emission, and strongly variable (the 05 stellar winds. This is the case for Wolf-Rayet stars, where the 
supergiant Cyg OB2 No. 9, for instance, was observed to vary correlation between non-thermal radio emission and binarity is 

firmly established, and binarity appears to be a prerequisite for 

Send offprint requests to: S. Van Loo, non-thermal emission (Dougherty & Williams 20001' ■ The situ- 

e-mail: svenvl@ast.leeds.ac.uk ation is less clear for O stars. Of the 16 currently known non- 



Hot stars are sources of both thermal and non-thermal radio 
emission. The thermal radiation is free-free emission from ther- 
mal electrons in the ionised stellar wind. As all hot stars have 
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Table 1. Binary status of non-thermal radio emitting O stars. 
The list is compiled from Bieging et al. J 19891 a), Persi et al. 
( fT985l b), Leitherer et al. ( IT9951 c), Benagha et al. dMHI d), 
Setia Gunawan et al. (EUU^ e), Benagha & Koribalski 
f) and Drake ( fT990l g). 



Name 


Spectral 


Binary 


NT 




Type 


Status 


Ref. 


HD 93 129 A 


02 If + ? 


BIN' 


f 


HD 93250 


03 V((f)) 




c 


9Sgr 


04 V 




a 


CD-47°4551 


05If 




d 


Cyg 0B2 No. 9 


05If 




a 


HD 150136 


05 Illn(f) + ? 


BIN^ 


d 


HD 168112 


05III 




a 


HD 15558 


05 III -H ? 


BIN^ 


a 


Cyg 0B2 No. 8A 


06 -H 05.5 


BIN'' 


a 


Cyg 0B2 No. 5 


06 f -H 07 f 


BIN^ 


b 


HD 124314 


06 V(n)((fO) + ? 


BIN*^ 


d 


Cyg OB2-335 


07 V 




e 


15 Mon 


07 Ve 


BIN** 


g 


HD 167971 


07.5 Vf-hO 


BIN' 


a 


5 0ri A 


09.5 II -H ? 


BIN' 


g 


a Ori AB 


09.5 V + ? 


BIN'" 


g 



' Walbom <2002k ^ Garmany et al. <I980> : ^ Garmany & Massey 
I198B ; " De Becker et al. O004l : ' Torres-Dodgen et al. (1991}; ^ 
Gies (MT); ^ Leitherer et al. (TW); ^ Gies et al. (T^; ' Harvey et 
al. {1987}; '" Bolton fT974l 



thermal radio emitting O stars, only 10 are confirmed binaries 
(TableQ. 

A priori, there is no need for non-thermal emitters to be 
binaries. The driving mechanism of hot-star winds is known 
to be subject to a line-deshadowing instability (Owocki & 
Rybicki 119841 Feldmeier 1200 It . Time-dependent hydrody- 
namical simulations including the line-deshadowing instability 
(e.g. Runacres & Owocki |2002> . show that it produces a wealth 
of shocks that in principle are capable of accelerating electrons 
to high speeds. We refer to these shocks as wind-embedded 
shocks. Non-thermal radio emission from hot stars was first in- 
terpreted as synchrotron emission from electrons accelerated 
in such wind-embedded shocks by White ( I1985> . Substantial 
improvements of this original idea were made by Chen Cl992i) 
and Chen & White ( fT^ . 

This is the third in a series of papers investigating the non- 
thermal radio emission from single stars. The general aim is 
to determine the properties that wind-embedded shock models 
must have to reproduce the observed non-thermal radio emis- 
sion. By this means, we are also able to shed some light on 
the fundamental question of binarity, as it will become clear to 
what extent the wind-embedded shock model is able to explain 
the observations. Such an indirect approach is necessary be- 
cause it is difficult, in general, to confirm observationally that a 
star is single. A failure of the wind-embedded shock model to 
explain the observations would strongly suggest that colliding 
winds are required to produce non-thermal radio emission. 



In Van Loo et al. J2004I hereafter Paper I), we intro- 
duced a simple parametrised model. The number distribution 
of synchrotron-emitting electrons was given by a power law, 
both as a function of momentum and distance. The emission 
region was assumed to extend continuously out to an outer ra- 
dius /?rnax- This assumption is open to criticism, because cool- 
ing mechanisms such as inverse Compton scattering can limit 
the emission to potentially thin layers behind the shocks (Chen 
I1992I I. This effect was included in Van Loo et al. J2005I here- 
after Paper II). In that paper we found that the emitting layers 
can indeed be quite narrow and are also few in number, be- 
cause the strongest shocks tend to dominate the emission. In the 
present paper, we further extend the model to take into account 
the fact that shocks gradually weaken as they move out with the 
flow. Specifically, we use results of time-dependent hydrody- 
namical simulations as input to our synchrotron-emission mod- 
els. 

The comparison of synchrotron spectra with hydrodynam- 
ical simulations of instability-generated structure has the addi- 
tional advantage of providing supplementary observational in- 
put for such simulations. Although the theoretical argument for 
the line-deshadowing instability is extremely strong (Owocki 
1991;), and the instability-generated structure is still the best 
explanation for, e.g., soft X-ray emission, black troughs, 
and clumping inferred from density-squared diagnostics (e.g. 
Owocki [l991> . the small-scale inhomogeneities it generates are 
not easy to detect. Therefore, any additional observational in- 
formation about such structure would be most welcome. 

The remainder of this paper is organised as follows. In 
Sect. 12 we introduce the model. In Sect. |3 we conclude that 
it fails to explain the Cyg OB2 No. 9 observations. We explore 
various "escape routes" in Sect. |3 Conclusions are given in 
Sect. 13 

2. The model 

2.7. Layered emission 

As the present model is based on the one from Paper II, 
we briefly review the properties of this latter model. The 
synchrotron emission is produced by electrons accelerated in 
wind-embedded shocks, produced by the line-deshadowing in- 
stability in a single star wind. The acceleration is not calculated 
ab initio, but a power law is assumed for the momentum distri- 
bution of electrons. 

In principle, the pressure from the relativistic particles can 
influence the structure of the shocks, and can even lead to their 
destruction (Drurv [T983> . Incorporating this influence in time- 
dependent hydrodynamical simulations is beyond the scope of 
the present work. We therefore neglect the effect of relativis- 
tic particles on the hydrodynamics (this is referred to as the 
"test-particle approach"). The result of incorporating this influ- 
ence in time-dependent hydrodynamical simulations would be 
to weaken the shocks, thereby making it more difficult to pro- 
duce the extended regions of non-thermal emission needed to 
explain the observations. As we shall find, it is already diffi- 
cult to produce the observed emission based on present shock 
strengths (Sect, \2.2l . so our conclusions will be all the more 
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valid if the feedback of relativistic particles on the hydrody- 
namics is included. 

We adopt the magnetic field expression of Weber & Davis 
J1967> . This means that, at large distances from the star, 
the magnetic field is almost perpendicular to the shock nor- 
mal. Jokipii (I1987> showed that this configuration (i.e. quasi- 
perpendicular shocks) is the most efficient for particle acceler- 
ation. Details are given in Van Loo (12005 > . 

The electrons are subject to a number of cooling mecha- 
nisms, most notably inverse Compton scattering and adiabatic 
expansion. Cooling has an important effect both at the shock 
front, and downstream. At the shock front, the maximum at- 
tainable energy is set by the balance between acceleration and 
cooling, imposing a high-momentum cut-ofi" on the distribu- 
tion. Downstream from the shock, the acceleration ceases and 
cooling rapidly removes the energy of relativistic electrons. 
The emission is then expected to come not from a continuous 
volume, but from a (possibly thin) layer behind each shock. 
Cooling also modifies the shape of the momentum distribution 
close to the cut-off. Away from the cut-off it retains its power- 
law shape. 

In Paper II, we also showed that the strongest shocks dom- 
inate the emission. This means that the number of synchrotron 
emission layers that contribute substantially to the observed 
emission is small. It is therefore natural to consider first the 
case of a single shock. A single-shock model is characterised 
by four free parameters: the position Rs of the shock, its veloc- 
ity jump Am and its compression ratio and the surface mag- 
netic field B,. The fluxes include a small but significant thermal 
contribution (25% of the total flux at 2 cm). For a complete dis- 
cussion of the model, we refer to Paper II. 

To describe the distribution behind the shock front, we need 
to know where the electrons were accelerated and how they are 
cooled. The relativistic electrons move away from the shock 
front with the gas. In terms of the velocity jump Am and the 
compression ratio the post-shock speed in the frame of the 
shock can be expressed as 



which follows directly from the Rankine-Hugoniot relations. 
As the gas flows a distance ARs away from the shock at a speed 
M2, the shock has travelled a much larger distance x AR^Uoa/ui, 
where Uca is the terminal speed. As Vao/ui ~ 100, the electron 
distribution in a fairly narrow layer behind an outer-wind shock 
reflects wind conditions over a wide range in distance. The 
outer electrons of a typical synchrotron-emitting layer (4 
wide, behind a shock at 500 R^), will have been accelerated at 
a distance ~ 100 from the central star. Of course, they are 
subjected to severe cooling during the time-interval needed to 
reach the outer wind. 

In Paper II, we assumed a constant shock velocity jump and 
compression ratio, mainly to keep the number of free parame- 
ters in the model as small as possible. However, hydrodynami- 
cal simulations show that both the shock velocity jump and the 
compression ratio decline with radius (see Fig. Q. The shock 
velocity jump appears as a strong scaling factor in this model. 



with the flux proportional to Am . The dependence of the syn- 
chrotron flux on the compression ratio x is somewhat more 
complicated, but is also important (Van Loo 120051 . Generally, 
the flux increases rapidly with x- In view of the wide range in 
radius reflected by the electron distribution, one can therefore 
expect that introducing the radial decline of the velocity jump 
and the compression ratio has important consequences for the 
models. 

Due to the free-free opacity of the dense stellar wind, any 
radio emission (thermal or non-thermal) that is emitted too 
close to the star will be absorbed. Roughly speaking, the emer- 
gent radio flux originates beyond the radius of unit free-free op- 
tical depth. This radius can be viewed as a kind of "radio pho- 
tosphere" of the stellar wind. Due to the wavelength-squared 
dependence of the free-free opacity, the size of the radio pho- 
tosphere increases with wavelength. For example, in the well- 
studied non-thermal radio emitting O star Cyg OB2 No. 9, the 
radio photosphere is 90 Rt at 2 cm (15 GHz) and 195 R, at 
6 cm (5 GHz). The emergent flux of a non-thermal emitter is 
determined by the tight interplay between non-thermal emis- 
sion and thermal absorption. This is what gives the models their 
great diagnostic value. A more detailed discussion of the radio 
formation region is given in Paper I. 

2.2. Hydrodynamical input 

Since the pioneering models of Owocki et al. J1988t . many 
hydrodynamical simulations of hot-star winds have included 
the line-deshadowing instability (e.g. Feldmeier il995i Dessart 
& Owocki 120021 Runacres & Owocki I2002t . These simula- 
tions show that within a few stellar radii of the surface, the 
flow becomes highly structured, with gas concentrated in dense 
clumps, and pervaded by strong shocks. These simulations are 
restricted to distances smaller than 100 which is not far 
enough to provide input for synchrotron-emission models. For 
simulations extending out to large radii, Runacres & Owocki 
J2005t introduced a periodic -box formalism for the outer evo- 
lution of structure in hot-star winds, which enabled them to 
follow structure out to distances of more than 1000 R,. 

In Fig.^we show the velocity jump and compression ratio 
of a typical strong shock in a periodic-box model as it moves 
out with the gas. The shock velocity jump steadily decreases 
from > 100 kms-' at 100 R, to * 55 kms"' at 250 R, and 
further to * 20 km s"' beyond 500 The compression ratio 
rises from ^ 2.6 at 100 Rt to 3.3 at 250 and then decreases 
to values » 2.2 beyond 500 Rt. The compression ratio increases 
between 100 and 200 Rt because the strong forward shock that 
we follow overtakes a few weak forward shocks and merges 
with them. 

As it takes some time to cool relativistic electrons below 
radio-emitting energies, the inner-wind accelerated electrons 
can contribute significantly to the radio emission. Furthermore, 
since the emissivity rapidly declines with radius due to the de- 
creasing shock velocity jump and compression ratio, the emer- 
gent flux is weighted in favour of electrons accelerated in the 
inner wind, and the radial dependence of the velocity jump and 
the compression ratio must be included in the model. For the 
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Fig. 1. Compression ratio (plus signs, upper panel) and ve- 
locity jump (plus signs, lower panel) of a typical shock in the 
periodic -box model of Runacres & Owocki J2005t . as a func- 
tion of radius. The solid lines represent the step function ap- 
proximations used in the code. 



sake of simplicity, we approximate the radial dependence of Am 
and^ by step functions of the form shown in Fig. [I] 



v(Hz) 

Fig. 2. Pure synchrotron flux for a single-shock model with 
declining Am and;^' given by the step functions in Fig.n(solid 
line), compared to a model with constant shock strength x - 
2.2 and Am = 18 kms"' (dashed line), i.e. a model where the 
history of the shock has not been taken into account. The shock 
is located at Rs = 600 R^. The surface magnetic field is - 
200 G. Stellar and wind parameters are as in Paper II. 



3. Results 

3.1. Single shock 

In a single-shock model, explicitly taking into account the ra- 
dial dependence of Am and x has an important eff'ect on the 
emergent flux. In Paper II we showed that a single shock lo- 
cated at ~ 600 Rt can explain the observed non-thermal emis- 
sion of Cyg OB2 No. 9, albeit with a much larger velocity jump 
(290 kms"') than predicted by hydrodynamical simulations. 
From Fig.[3we derive that a typical shock at that distance has a 
compression ratio - 2.2 and a velocity jump Am = 18 kms"'. 
In Paper II we assumed that the shock strength is the same ev- 
erywhere in the wind, whereas the hydrodynamical simulations 
predict a weakening of the shock. It is to be expected, for a 
weak outer-wind shock, that its past history as a stronger inner- 
wind shock affects the energy distribution of the accelerated 
electrons behind it. 

In Fig. |2] we compare the pure synchrotron flux (without 
the thermal contribution) generated by such a single shock as- 
suming a constant shock strength throughout the wind (dashed 
Une), and the flux emitted by the same shock when its history 
is taken into account (solid line). This last model produces sig- 
nificantly more emission (by a factor of ~ 200). This is be- 
cause, when the shock was closer to the star and consequently 
stronger, it produced a population of electrons that are still suf- 
ficiently energetic, in spite of strong inverse-Compton cooling, 
to produce synchrotron emission. 

When the history of the shock is taken into account, a single 
shock can still explain the observed non-thermal spectrum of 
Cyg OB2 No. 9. To fit the observed flux, we need to adjust the 



magnitude of the shock velocity jump. We do this by defining 
the velocity jump at radius r as 

AM(r) = Au^^^f(r), (2) 

where f{r) is the velocity jump from Fig. [2 divided by its high- 
est value (i.e. f(r) < 1), and AM,nax is a scaling factor that sets 
the magnitude of the velocity jump. 

In Fig. 13 we show a single-shock model that fits the ob- 
served spectrum. For simplicity, the compression ratio x is 
taken to be 4 at all radii. The value of AM^ax is 1050 kms"', 
corresponding to a local velocity jump at Rs - 602 R, of 
164 km s"' . This is lower than the value determined in Paper II 
(due to the significant amount of inner wind accelerated elec- 
trons in the emission layer) but still an order of magnitude 
higher than predicted by hydrodynamical simulations. With a 
decreasing x, the value obtained for Am would be even higher. 
The only way to fit the observations with smaller velocity 
jumps is to consider multiple shocks. 

3.2. Multiple shocks 
3.2.1 . Constant^ and Am 

The synchrotron emission in our model is proportional to the 
third power of the shock velocity jump (Paper II). This allows 
us to estimate the flux of a multiple-shock model with a con- 
stant compression ratio and velocity jump from a single-shock 
model. The condition for a collection of identical shocks 
with velocity jump Am, distributed throughout the wind, to pro- 
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Fig. 3. A single-shock model (solid line) that fits the observa- 
tions (diamonds) of Cyg OB2 No. 9. The vertical stripes are the 
error bars. The shock velocity jump is IS.u{r) - Au^-^s^fir), as 
defined in Eq. (|2ji. The fit parameters are Awmax - 1050 km s"', 
i?s = 602 ;if = 4 and = 200 G. 

duce the same flux as a single shock with velocity jump AMsingie 
is simply expressed as 

NAu' = Am^.^„i,. (3) 

This offers a way of lowering the velocity jump needed to ex- 
plain the observations: one just adds shocks until the veloc- 
ity jump is in the range of the hydrodynamical predictions. 
In Paper II, we fitted the observations by a single shock with 
Ausingie = 290 kms"'. Using the above equation we estimated 
that 200 shocks with Am = 50 kms ' should produce the same 
flux. 

When the shock velocity jump and/or compression ratio 
decrease as a function of radius, as predicted by hydrody- 
namical simulations, Eq. Q can no longer be applied. In that 
case, we find that our model cannot reproduce the observations 
(Sect. l3T2!2t . To understand why this is so, it is useful to con- 
sider first the cumulative flux for a model with a constant shock 
velocity jump and compression ratio. (The cumulative flux at 
a position r in the wind is defined as the emergent flux due 
to all the shocks between the stellar surface and r.) In Fig. |3 
we show the cumulative fluxes at different radio wavelengths 
for an equidistant distribution of shocks, with a shock at every 
12 Rt (The influence of the adopted shock spacing is discussed 
in Sect. 14. U . All shocks — 3-7 and Am - 71 kms"'. The 

surface magnetic field is B, = 200 G. Only the synchrotron 
flux is shown. The "observed" synchrotron flux is obtained by 
subtracting the theoretical thermal flux (using the Wright & 
Barlow formalism) from the observed radio flux. The cumu- 
lative synchrotron flux at 2 cm (solid line) is zero for distances 
below ~ 50 7?,, where all synchrotron emission is absorbed 
by the free-free opacity of the wind, and then rises sharply. The 
flux at 6 cm (dashed) begins to rise further out, due to the larger 
free-free optical depth. For the same reason, the 20 cm begins 
to rise even later. 
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Fig. 4. Cumulative fluxes for a model with constant compres- 
sion ratio and velocity jump, at 2 cm (solid line), 6 cm (dashed) 
and 20 cm (dashed-dotted). The diamonds are the observed 
synchrotron fluxes and also the error bar on the observations 
is given. The model parameters are Am = 71 kms"', x - ^-1 
and B, = 200 G (Am and^ are constant in the wind). A shock 
is present every 12 



The cumulative flux increases more steeply at larger wave- 
lengths. The reason is twofold. First, there is the intrinsic spec- 
trum of synchrotron emission: the flux at larger wavelengths is 
produced by less energetic electrons, of which there are more 
by virtue of the shape of the momentum distribution. Second, 
inverse-Compton cooling is less severe at larger wavelengths, 
so that the synchrotron-emitting layer is wider. 

3.2.2. Declining and Am 

Let us now see what happens when the shock velocity jump 
and compression ratio decrease as a function of radius, as pre- 
dicted by hydrodynamical simulations. We can no longer use 
Eq. (|3} to determine how many shocks are needed, but plotting 
the cumulative flux allows us to judge when we have a good fit 
to the observations, without requiring the number of shocks to 
be pre-specified. 

In Fig.|5]we show the cumulative flux for a model that is 
similar to the model shown in Fig.|?] except for the decreasing 
shock strength. The compression ratio of the shocks is given 
by Fig.[2(top panel). The velocity jump is given by the lower 
panel of the same figure, but a scaling factor AMmax is intro- 
duced, as in Sect. 13.11 to fit the observed flux at 6 cm. In the 
model shown, AM^ax = 194 kms"', which represents a 70% 
increase with respect to the hydrodynamical predictions. The 
cumulative fluxes are zero below the radio photosphere, and 
then rise steeply. Contrary to the model with constant velocity 
jump and compression ratio, there is no discernible difference 
in the rate of increase of the fluxes at 2 and 6 cm. If the shock 
strength were constant, the 6 cm flux would rise faster than the 
2 cm flux (see Sect. 13 .2.11 . This is now compensated by the fact 
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Fig. 5. Cumulative fluxes for a model with decreasing com- 
pression ratio and velocity jump, at 2 cm (solid line), 6 cm 
(dashed) and 20 cm (dash-dotted). The diamonds are the ob- 
served synchrotron fluxes and also the error bar on the ob- 
servations is given. The model parameters are B, - 200 G, 
Awmax = 194 kms"' and;if is taken from Fig.[n A shock is 
present every 12 

that the 2 cm flux is formed at smaller radii, where the shocks 
are stronger. 

Both the 2 and 6 cm fluxes level off at large radii, because 
the compression ratio and velocity jump become so small that 
the emitted synchrotron radiation is negligible. The fluxes at 2 
and 6 cm reach approximately the same asymptotic value, re- 
sulting in a flat synchrotron spectrum. After adding the thermal 
contribution, the emergent spectrum has a positive slope, con- 
trary to observations. At 20 cm, the model significantly under- 
estimates the observation. Extending the model to larger radii 
would increase the 20 cm flux somewhat, but would have no 
effect on the 2 and 6 cm flux. 

In summary, we find that a model where the velocity jump 
and compression ratio are taken from hydrodynamical simula- 
tions does not reproduce the observed radio spectrum, because 
the strong shocks close to the star contribute more to the flux 
at shorter wavelengths (where the free-free opacity is smaller) 
than at longer wavelengths. Any model where the emissivity 
is a strongly decreasing function of radius will therefore have a 
tendency toward a positive spectral index. Conversely, to repro- 
duce the observed negative spectral index, any model should 
have an emissivity that is a weakly decreasing function of ra- 
dius. 

4. Discussion 

The hydrodynamical simulations from which we take the shock 
velocity jump and compression ratio are subject to considerable 
uncertainty. Foremost, they are performed in one dimension, 
whereas the stellar wind is expected to contain a fair amount 
of lateral structure (e.g. due to Rayleigh-Taylor instabilities, 
see Dessart & Owocki |2003> . Moreover, no attempt has been 



made to model a particular star. The results of these simulations 
should therefore not be taken at face value. To further quan- 
tify the discrepancy between the hydrodynamical predictions 
and the radio observations, let us see by which modifications 
the two may be reconciled or, in other words, which "escape 
routes" may lead us away from the discrepancy, and whether 
any of these escape routes is physically plausible. 

4.1. Shock spacing and outer boundary of the 
synchrotron-emitting region 

In Sect. |3lwe adopted a shock spacing of 12 This value is 
somewhat arbitrary but is based on the spacing found in hydro- 
dynamical simulations, and on the fact that only the strongest 
shocks will contribute substantially to the synchrotron emis- 
sion. Reducing the spacing from 12 /?« to 4 Rt, which is about 
the smallest distance between two adjacent outer wind shocks 
in the hydrodynamical simulations, does not change the slope 
of the emergent synchrotron spectrum. The only effect of the 
closer spacing is a lower value of Au^^x, due to the larger num- 
ber of shocks in the wind. 

To produce the observed negative spectral index, the radial 
decrease of the synchrotron emission should be weak. Such a 
slow decrease can be obtained by a closer shock spacing in the 
outer wind than in the inner wind (to counteract the decrease 
of the shock strength). Physically this corresponds to a slow- 
ing down of the shocks as they move out with the wind. Such 
a slowing down can indeed be expected from momentum con- 
servation. However, the associated change in velocity is mini- 
mal. Both the shell and the material it sweeps up are moving at 
roughly the terminal speed, with velocity differences of the or- 
der of a few times the sound speed. We find that, to reproduce 
the observed spectral index, the spacing of the shocks in the 20 
cm formation region should be less than half of that in the 2 
cm formation region. This corresponds to a decrease in shock 
speed (measured in the stellar rest frame) of a factor of two or 
more, which is not predicted either by detailed hydrodynam- 
ical simulations or by the momentum conservation argument 
above. 

The adopted outer boundary of the synchrotron-emitting re- 
gion could in principle influence the spectral index predicted by 
the model. However, the fact that the cumulative flux in Fig.|5] 
levels off at large radii shows that the spectral shape cannot be 
changed by extending the synchrotron-emitting region to larger 
radii (as the contribution at large radii becomes negligible). 

We conclude that our model cannot be made to produce a 
negative spectral index by increasing the number of shocks in 
the wind, either by changing the shock spacing or by adopting 
a larger outer boundary of the synchrotron-emitting region. 

4.2. Rate of decline of x ^nd Au 

To reproduce the observed spectra, any model should have an 
emissivity that does not decline too rapidly as a function of ra- 
dius. As the emissivity is directly related to the compression 
ratio and velocity jump of the shocks producing the emission, 
making these properties a weaker function of radius is a nat- 
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Fig. 6. Cumulative fluxes for a model with a "flattened" Au{r), 
at 2 cm (solid line), 6 cm (dashed) and 20 cm (dash-dotted). 
The diamonds are the observed synchrotron fluxes and also 
the error bar on the observations is given. The ratio between 
the maximum and minimum Am is 3. The other parameters are 
Awmax = 835kms"',;^f = 2 andB* - 200G. Ashockis present 
every 12 

ural way of obtaining a more slowly declining emissivity. As 
it is difficult to disentangle the effects of a declining velocity 
jump and a declining compression ratio, we adopt a constant 
compression ratio and see how weak a function of radius the 
velocity jump should be to produce a negative spectral index. 
Therefore, we make AM(r) a weaker function of r, by replacing 
/(r) in Eq. Q with fir), defined as 

/>)= l-e[l-/(r)] 0<e<l (4) 

where e - I corresponds to the original /(r) and e = to a 
constant velocity jump. In Fig.|6lwe show the cumulative flux 
for a model where AM(r) has been "flattened" in this manner 
and a constant compression rdX\o x - 2. The ratio between the 
maximum and minimum values of Am is 3 (corresponding to 
e = 0.79), as opposed to 6.4 in Fig.^ This model explains the 
observations at all wavelengths, although the required Am^^x is 
much higher than the hydrodynamical predictions. The obser- 
vations can also be reproduced with a higher compression ratio 
X = ^ and an even flatter AM(r), with a ratio of 1.5 between its 
maximum and minimum. 

In these models, we assumed a constant compression ra- 
tio, whereas the hydrodynamical simulations predict a decline 
(Fig. With a declining compression ratio, AM(r) must be 
even flatter. 

Indeed, we find that, with;\f(r) as defined in Fig.^ Am must 
be constant to fit the observations. If x{f) is flattened in the 
same way as was done above for AM(r) (Eq. Q), a fit can be 
obtained if the ratio between maximum and minimum is 1.3 or 
less for both Am andx- 

The upshot of the above is that the radio spectrum of 
Cyg OB2 No. 9 can only be explained with our model with 



a modest (< 30%) decrease of the velocity jump and compres- 
sion ratio. The range in radius contributing to the emergent flux 
at 2, 6 and 20 cm is ~ 30 - 1 000 Although not absolutely 
impossible, it is unlikely that shocks would have such a small 
decrease in velocity jump and compression ratio over such a 
large distance. 

4.3. Re-acceleration 

Our results show that it is difficult for a wind-embedded shock 
model to reproduce an observed non-thermal radio spectrum, 
if the radial decline of the shock strength is in agreement with 
hydrodynamical predictions. This may seem at odds with some 
previous models that successfully reproduce the observed non- 
thermal spectra. The most advanced of these models is given in 
Chen iTWI\ and Chen & White iT^^ . Their model can suc- 
cessfully reproduce the observed non-thermal spectrum of both 
9 Sgr and Cyg OB2 No. 9. It includes adiabatic and inverse- 
Compton cooling. As a result, electrons cannot travel very far 
from the acceleration site, and the electrons have to be acceler- 
ated in situ in the radio-emitting region. The emission is then 
also confined to layers behind the shocks. In that respect, the 
Chen & White ( II994> model is very similar to the model pre- 
sented in the present paper. 

The fundamental difference between Chen & White and the 
present model is re-acceleration of electrons. If a fluid element 
goes through a number of successive shocks as it moves out- 
ward, then relativistic electrons in the outer wind will, on av- 
erage, have passed through more shocks than inner-wind elec- 
trons. Consequently, these outer- wind electrons are more en- 
ergetic than they would have been at the same position in the 
wind, had they passed through only one shock. Therefore, re- 
acceleration substantially slows down ffie outward decrease of 
the synchrotron emissivity and is expected to facilitate a neg- 
ative spectral index. Indeed, with a simplified model of re- 
acceleration, where we artificially slow down the radial de- 
crease of the relativistic electron density, we easily reproduce 
the observations. 

Whether or not re-acceleration occurs depends on the 
adopted hydrodynamical model. In Chen & White ( II994> . the 
hydrodynamical model is from Lucy JI982> . supplemented 
with a snowplough model to describe structure in the outer 
wind (Chen I1992> . In Lucy (I1982> . all shocks are forward 
shocks. As was pointed out by White (White & Chen lI994l 
discussion), a fluid element in the Lucy model is accelerated 
to the terminal speed only in the shock jumps. Therefore, any 
given fluid element goes through a large number of shocks. 
Relativistic electrons are then accelerated at one shock front, 
cool as they travel downstream to the next shock, and are re- 
accelerated there. A population of new, freshly accelerated rel- 
ativistic electrons is also created at each shock. 

Lucy's model is entirely phenomenological, however. 
When the line-deshadowing instability is introduced in a model 
that solves the time-dependent equations of hydrodynamics, 
the physical picture of wind structure is quite different. The 
fundamental consequence of the instability is the omnipresence 
of strong rarefaction waves of fast, rarefied gas. These rarefac- 
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tions are terminated by a reverse shock (much as when gas runs 
into a solid wall), forming the inner boundary of a dense shell. 
As this dense shell is impacted by faster shells, it develops a 
forward shock on the outside. The result is that shocks occur 
in pairs at either side of a dense shell. Gas is fed through the 
shocks into the shell and then remains trapped there. A fluid 
element therefore only encounters a single shock, making re- 
acceleration impossible. 

The simulations on which the above analysis is based are 
one-dimensional. In a three-dimensional situation the trapping 
of the gas is less absolute, as gas can in principle escape side- 
ways from a shell fragment or blob and encounter another 
shock. It is hard to make a quantitative estimate of the impor- 
tance of such an effect. As the radiative driving and the associ- 
ated instability are strongest in the radial direction, one can ex- 
pect the azimuthal velocity to be relatively small, even though 
instabilities such as Rayleigh-Taylor are expected to disrupt the 
lateral coherence of the structure. One situation where the az- 
imuthal velocity could become more important is the case of 
rarefied gas flowing around a dense blob. However, this gas has 
a low density and should therefore not contribute much to the 
emission. It therefore seems likely that most of the dense gas 
will see few shocks, which impairs the main mechanism by 
which Chen & White J1994> were able to produce a negative 
spectral index. 

4.4. Smaller radio photosphere 

The main impediment for a strongly decreasing synchrotron 
emissivity to produce a negative spectral index is the differ- 
ence in the radio photosphere at 2 and 6 cm, which causes 
the 2 cm flux to benefit from strong emission that is inacces- 
sible to the 6 cm flux. This difference in radio photosphere in- 
creases with density, and hence with the mass-loss rate. It can 
therefore be expected that a lower mass loss rate may produce 
a negative spectral index even with a strongly declining syn- 
chrotron emissivity. In Fig. we show the cumulative flux for 
M = 4.0x10"'' M0yr-', which is a factor 5 lower than the ob- 
served Ha mass-loss rate we adopted in the rest of this paper. 
By restricting the outer boundary of the synchrotron-emitting 
region to 430 we find that the model reproduces the ob- 
served radio spectrum. 

There is no reason to assume, however, that the mass- 
loss rate of Cyg OB2 No. 9 can simply be lowered by a fac- 
tor of five. Although the mass-loss rate is extremely high for 
an O star, Cyg OB2 No. 9 is also extremely bright (L, = 
1 .74 X 10^ L0), and it obeys the normal wind-momentum lumi- 
nosity relation for O stars. It is widely accepted that observed 
O-star mass-loss rates might be overestimated by a factor of a 
few due to clumping (e.g. Massa et al. 120031 Hillier et al. l2003> . 
The hydrodynamical simulations on which we have based our 
model obviously predict a large amount of clumping. However, 
it is important to realise that a reduced mass-loss rate due to op- 
tically thin clumps does not reduce the free-free optical depth 
if the clumped model is to produce the same thermal flux (see 
Appendix IXt. In other words, assuming a clumped wind does 
not provide a way to reconcile the hydrodynamical predictions 
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Fig. 7. Cumulative fluxes for a synchrotron model with an ar- 
tificially low mass-loss rate, at 2 cm (solid line), 6 cm (dashed) 
and 20 cm (dash-dotted). The diamonds are the observed syn- 
chrotron fluxes and also the error bar on the observations is 
given. The model parameters are M = 4.0 x 10"^ Mgyr"', 
Awmax = 152 kms"',;if = 3 and = 200 G. We extrapolated 
Am inward, so that, for r < 30 7?*, Am = 1.74 AMn,ax and, for 
30 7?* < r < 100 Rt, Am = 1 .39 AM^ax- A shock is present every 
MR,. 

with the radio observations. Other mechanisms such as poros- 
ity (e.g. Owocki et al. I2004> . where wind material is packed 
together also in the lateral direction, may play a role when in- 
dividual clumps become optically thick. Owocki et al. (I2004> 
show that for optically thick clumps the opacity is effectively 
reduced. Porosity is characterised by an additional parameter 
(the clump size or equivalent) which complicates the relation 
between the clump properties and the radio flux. It is as yet un- 
clear to what extent hot-star winds are affected by porosity, if 
at all. 

The 04 main-sequence star 9 Sgr has a much lower mass- 
loss rate (M = 2.4 x 10"^ Moyr"') flian Cyg OB2 No. 
9. Therefore, one could expect its spectrum to be compati- 
ble with the predicted decline of the shock strength. Figure |8l 
shows the cumulative flux at 2, 6 and 20 cm for a synchrotron 
model of 9 Sgr. The adopted parameters are T^s - 43 000 K, 
R, = 16 Rq, d = 1.58 kpc, = 2 950 kms"', M = 
2.4 X lO"*" Mq yr"', and L, = 8.1 x 10^ Lq. The velocity jump 
and compression ratio were taken from Fig.^ The model does 
not fit the observations, because the theoretical synchrotron 
spectrum between 6 and 20 cm is essentially flat, contrary to 
the observed negative spectral index. In this sense the same 
problems as for Cyg OB2 No. 9 appear, but now at larger wave- 
lengths. 

5. Conclusions 

We have included the radial decline of the shock strength, 
as predicted by hydrodynamical models, in our synchrotron- 
emission model for single O stars, and applied it to the widely 
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presumed single O star Cyg OB2 No. 9, which is a synchrotron 
emitter. We find that the decline of the shock strength no longer 
allows a wind-embedded shock model to reproduce the ob- 
served negative spectral index of non-thermal emission. We 
have investigated a number of "escape routes" by which the 
hydrodynamical predictions might be reconciled with the ra- 
dio observations. We find that a negative spectral index can 
be produced by a slower decline of the compression ratio and 
the velocity jump, by the re-acceleration of electrons in many 
shocks or by adopting a lower mass-loss rate. However, none 
of these escape routes are physically very plausible. In particu- 
lar, re-acceleration by feeding an electron distribution through 
a number of shocks, is in contradiction with modern hydrody- 
namical calculations. This leads us to the conclusion that there 
is probably no viable way to produce the observed spectrum 
based on hydrodynamical simulations. 

Our model obviously has its limitations. Most important, 
perhaps, is the use of one-dimensional (ID) hydrodynamical 
simulations. Two-dimensional simulations have been presented 
for instability-generated structure in the inner wind (Dessart & 
Owocki^OOS ). Overall, these models are characterised by a re- 
duced dumpiness of the wind and a larger (radial) velocity dis- 
persion. It is not feasible at present to perform 2D-simulations 
of the wind out to the distances required for synchrotron- 
emission models. It cannot be excluded that such models might 
show a slower decrease of the shock strength. It is to be ex- 
pected though that as radiative driving becomes negligible at 
distances of a few tens of stellar radii (Runacres & Owocki 
I2002t . both the compression ratio and the velocity dispersion 
should eventually decrease even in a 2D model. Furthermore, 



2D and 3D hydrodynamics allows for many additional insta- 
bilities (such as Rayleigh-Taylor) to disrupt the structure. It is 
impossible at present to make even an educated guess what a 
structure of such complexity might mean for the efficiency of 
synchrotron emission. 

Based on the current hydrodynamical simulations, which 
admittedly are open to question in their details, we suspect that 
the observed non-thermal radio emission from O stars cannot 
be explained by wind-embedded shocks associated with the 
instability of the line-driving mechanism. The most plausible 
alternative to emission by wind-embedded shocks is emission 
from shocks associated with colliding stellar winds. That would 
imply that all non-thermal O stars should be members of a bi- 
nary or multiple system. 

This conclusion is supported by the recent finding (Van Loo 



120051 Blomme et al. 2005, Rauw et al. 120051 . based on both 
period analyses of radio light curves and spectroscopic evi- 
dence, that three presumably single stars with non-thermal ra- 
dio emission are probably binaries. Continued systematic radio 
and spectroscopic monitoring may lead to the same conclusion 
for the remaining non-thermal radio emitters that are presumed 
to be single O stars. 
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Appendix A: Effect of clumping on the free-free 
optical depth 

The free-free optical depth is usually calculated assuming a 
smooth wind (Wright & Barlow II975> and the effect of den- 
sity inhomogeneities (clumping) is neglected. However, the 
line-deshadowing instability packs the stellar wind material to- 
gether in dense clumps (e.g. Runacres & Owocki 2002). In 
principle, because the free-free optical depth is proportional to 
the density squared, one expects clumping to have an influence 
on the radio photosphere of the stellar wind. 

This is indeed the case when we take a given amount of 
stellar wind material and consider the effect of clumping it. 
Due to the enhanced opacity of the clumped wind, a given op- 
tical depth (integrated from a location in the wind out to the 
observer) is reached closer to the observer, and as a result the 
radio photosphere of the stellar wind is larger. 

However, if we compare a smooth wind and a clumped 
wind that produce the same radio flux, the effect of the en- 
hanced optical depth is cancelled exactly by the fact that the 
clumped wind requires a lower mass loss rate to explain the 
observations. If we assume a uniform, random distribution of 
optically thin clumps, then we can use the description of clump- 
ing proposed by Abbott et al. JI98I> . From their Eq. (9) we 
see that the optical depth has the following dependence on the 
clumping factor /d and the mass-loss rate 

Ty^M^U- (A.1) 
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(A succinct explanation of the clumping factor is given in 
Runacres & Owocki J2002l i). The emergent radio flux scales 



Fy OC 



(A.2) 



(Abbott et al. ll981l their Eq. (10)). Therefore, a clumped wind 
and a smooth wind that produce the same radio flux will also 
have the same optical depth. Consequently, a given optical 
depth is reached at the same radius. This means that, in the 
limit of optically thin clumps, clumping has no effect on the 
values we derive for the radio photosphere of the stellar wind. 

When individual clumps become optically thick, the Abbott 
et al. formalism is no longer valid. In that case, porosity ef- 
fects can play a role. Owocki et al. 12004} show that for op- 
tically thick blobs the opacity is effectively reduced. Porosity 
is characterised by an additional parameter (the clump size or 
equivalent) which complicates the relation between the clump 
properties and the radio flux. 
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